Electrostriction and guidance of sound by light in optical fibers 
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We investigate the generation of phonon wavepackets in optical fibers via electrostriction from coherent 
optical waves. Solving the elastodynamic equation subject to the electrostrictive force, we are able to reproduce 
experimental spectra found in standard and photonic crystal fibers. We discuss the two important practical 
cases of forward interaction, dominated by elastic resonances of the fiber, and backward interaction, for which 
an efficient mechanism of phonon guidance by light is found. The latter result describes the formation of the 
coherent phonon wavepacket involved in stimulated Brillouin scattering. 



Electrostriction (ES) describes the coherent generation of 
acoustic phonons from the interference of two frequency- 
detuned optical waves. The effect is closely related to Bril- 
louin light scattering (BLS), or scattering of an incident pho- 
ton by an acoustic phonon of a solid material, accompanied 
by a frequency shift of the scattered photon. The combina- 
tion of BLS and ES explains stimulated Brillouin scattering 
(SBS), an efficient three-wave interaction involving two op- 
tical waves and a phonon wavepacket satisfying both energy 
and momentum conservation. SBS is a fundamental limit- 
ing phenomenon for signal transport over optical fibers [1]. 
Interactions of photons and acoustic resonances of photonic 
crystal fibers (PCF) [2] and optical waveguides [3] have re- 
ceived much attention, for instance in view of obtaining strong 
optoacoustic interactions [4] or of reducing noise in quan- 
tum optics experiments [5]. ES in optical fibers is not eas- 
ily observed, since while the incident and the scattered optical 
waves are directly accessible to experiment, the intervening 
acoustic phonons are only witnessed via their induced effects 
on light. As a result, most theories eliminate phonon dynamics 
to concentrate on the time evolution of the amplitudes of opti- 
cal waves. This process, however, can only be performed un- 
der strong assumptions on the nature of the phonon wavepack- 
ets generated by ES. Guided acoustic wave Brillouin scat- 
tering (GAWBS) theories, for instance, assume the acoustic 
modal shapes and resonant frequencies to be known before- 
hand, so that scattering cross-sections can be computed from 
overlap integrals [6-8]. Modal approaches, however, give no 
insight into the mechanism of selection of particular phonons 
from the electrostrictive force. 

In this Letter, we propose a direct calculation of phonon 
wavepackets induced by electrostriction in optical fibers, bas- 
ing only on knowledge of the fiber geometry and of the par- 
ticular incident guided optical waves. Specifically, an elasto- 
dynamic equation is solved for a fiber with an arbitrary cross- 
section subject to an optical force. By comparing with ex- 
perimental spectra obtained with standard optical fiber and 
photonic crystal fiber with a solid core, we show that the 
computed phonon wavepackets as a function of detuning fre- 
quency explain both backward and forward ES and provide an 
estimation of the efficiency of the process. Under the phase- 
matching conditions of SBS, furthermore, our results imply 
that the phonon wavepacket generated by electrostriction is 
actually guided by light. 



Standard SBS models for optical fibers are based on a rather 
strong plane-wave approximation [9, 10] and on the picture 
of ES-generated pressure or density fluctuation waves [1]. In 
contrast to this simple picture, theory and experiment have 
shown the relevance of including both shear and longitudinal 
elastic waves together with a precise description of the struc- 
ture [2, 6, 7, 11]. Let us consider an optical fiber with an 
arbitrary cross-section. We assume the total incident optical 
field results from the superposition of two frequency-detuned 
guided modes 

E(r,z;t) = EW^y^it-^z) + E (2)( r y(u 2 t-k 2Z )^ (1) 

with angular frequencies loi(2) and axial wavevectors fci(2). 
This optical distribution induces an optical force given 
by divergence of the symmetric ES stress tensor T?? = 
—eoXkiijEkEi , with the rank-4 susceptibility tensor XkUj = 
£k m £inPmnij and p mnij the elastooptic tensor. e is the per- 
mitivity of vacuum. The force term with detuning frequency 
lj = lo\ — L02 is proportional to E^'E^** exp(z(wi — kz)) 
with k — ki — &2- Because we are considering the generation 
of elastic waves with low frequency with respect to optical 
frequencies, lu <C uji.2- If the two optical waves are propagat- 
ing in the same direction, a situation typical of GAWBS [6], 
k and we speak of forward ES. If they are propagating in 
opposite directions, which is the SBS situation, k w 2ki and 
we speak of backward ES. Because of the particular source 
term, the following ansatz is assumed for the displacements 
of the ES-generated elastic wave 

u i {r,z;t)=u i (r)e i ^ t - kz \ (2) 

where the transverse dependence Ui{r) is the unknown of the 
model. The elastodynamic equation 

P~Q^T ~ ( C ijklU k j), 3 = -Tifj, (3) 

with Cijki the rank-4 tensor of elastic constants, is written in 
variational form by left-multiplying by virtual displacement 
field Vi and integrating over the cross-section, S. Further ap- 
plying Green's theorem, we get 

- cj 2 / pv*u t + / v-jdjkiukj = / drv*jT?j, (4) 
Js Js Js 
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which amounts to the theorem of virtual work for the ES 
force. For practical computations, the Galerkin nodal finite 
element method (FEM) is then employed to transform the in- 
tegral equation into the linear system 

(K(k) - uj 2 M)U = X(k)T es , (5) 

with mass matrix M, stiffness matrix K(k) = K a + kKi + 
k 2 K 2 , and X(k) = X + kX x . U is the vector of nodal dis- 
placements ui and T es is the vector of electrostriction stress 
tensor value at nodal points, = —eaXkUjEf, ' E\ 

Equation (5) is our central result. Its solution as a function of 
frequency detuning gives the rigorous distribution of displace- 
ments within the waveguide cross-section and can be directly 
compared to experiment, as we show next. The ES stress 
tensor is uniquely defined by the optical modal distribution. 
Matrices K (k) and X(k) are polynomials in the wavevector 
detuning k, thus unifying in the same formula forward and 
backward ES, as well as any intermediate situation. 

As a test for our model, we selected a 1 km long single 
mode fiber (SMF28) and a 400 m long large core PCF with 
similar optical propagation properties. Single mode guidance 
is achieved by total internal reflection: SMF28 has a Ge-doped 
core with larger refractive index than surrounding undoped 
silica, while PCF has a pure-silica core surrounded by a holey 
structure with lower effective index. The silica PCF shown 
as an inset in Fig. 2(a) has a hole diameter of 4.6 /im and 
an air filling ratio d/A= 0.58, resulting for the fundamental 
mode in effective index n e g = 1.441 and optical effective 
area A eS = 70 /im 2 . The SMF28 fiber has n cS = 1.446 and 
A e g = 78.3 /im 2 , while the core diameter is 8.2 /im. 

ES-generated phonons inside an optical fiber are hardly 
observable directly. They can, however, be observed indi- 
rectly in a manner similar to spontaneous Brillouin gain spec- 
trum measurement or Brillouin sensing, by monitoring light 
diffracted from a coherent pump wave. A heterodyne detec- 
tion experimental setup was used for investigating backward 
ES [12]. To observe forward ES, we used a fiber loop mir- 
ror interferometer [8]. The measured spectra for the SMF28 
fiber are show in Fig. l(a,b). The forward spectrum has many 
sharp peaks that gradually vanish for frequencies above 600 
MHz. The peaks are known to be directly related to elastic res- 
onances of the fiber. The backward spectrum shows a single 
Lorentzian peak with a linewidth of 27 MHz and a central fre- 
quency 10.85 GHz. The PCF case shown in Fig. 2(a,b) looks 
similar, though resonance frequencies are slightly different. 
Forward ES resonances in this case are known to be related to 
the holey structure rather than to the external cladding of the 
fiber. 

We then proceeded to model experimental results using 
Eq. (5). A FEM mesh covering the full fiber cross-section was 
first created, with the outer boundaries left free. The mesh of 
the PCF was extracted from the scanning electron microscope 
image shown as an inset in Fig. 2(b). Fundamental optical 
modes are first obtained and normalized for unit transported 
power. ES stress distributions are displayed in Figs. 1(e) and 
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FIG. 1. Comparison between experiment and ES model in single 
mode fiber (SMF28). (a) Experimental and (c) numerical forward 
spectra (k = 0). (b) Experimental and (d) numerical backward spec- 
tra (k — 2ki). (e) Spatial distribution of the ES force generated by 
the fundamental optical mode at a wavelength of 1550 nm with 1 W 
optical power. The black circle limits the germanium doping core 
area. (f,g) Computed kinetic energy density of phonon wavepackets 
at 178 MHz (forward ES) and 10.85 GHz (backward ES). 



TABLE I. Independent material constants for isotropic silica [13]. 
Values for 0.36%-GeO2 -doped silica given in parenthesis are de- 
duced from [14]. 

Elast. const., GPa c u = 78 (76) c 12 = 16 (16.15) c 44 = 31 (29.9) 
Photoelast. const. p n = 0.12 p 12 = 0.27 p 44 = -0.073 
Mat. dens., kg/m 3 p = 2203 (2254) 
Refractive index n = 1.444 (1.4492) at A = 1550 nm 



2(e). For both fibers, the ES stress is confined in the core 
roughly as the square of the optical field. Next, the ES-driven 
elastodynamic equation was solved for the displacement of 
the elastic wave by imposing the phonon wave vector k and 
scanning the detuning frequency w. Independent material con- 
stants considered for silica are shown in Table I. To include 
the phonon lifetime, elastic losses were incorporated in the ES 
model by considering a complex elastic tensor whose imagi- 
nary part is a constant viscosity tensor times frequency [15]. 
This loss model is compatible with the usual assumption that 
the Qf product is a constant for a given material, with Q the 
quality factor and / the frequency. This is the only adjustable 
parameter in the model. The value Qf — 5.10 12 Hz was se- 
lected to fit the SBS linewidth of SMF28 at 1 1.01 GHz. 

Calculated forward (k — 0) and backward (A: = 2k{) ES 
spectra are shown in Fig. 1 and Fig. 2 for SMF28 and PCF, 
respectively. These spectra are obtained by evaluating the ki- 
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FIG. 2. Comparison between experiment and ES model in large core 
PCF. (a) Experimental and (c) numerical forward spectra (k — 0). 
The inset in (a) shows a cross-section of the fiber, (b) Experimental 
and (d) numerical backward spectra (k — 2ki). (e) Spatial distribu- 
tion of the ES force generated by the fundamental optical mode at 
a wavelength of 1550 nm with 1 W optical power. (f,g) Computed 
kinetic energy density of phonon wavepackets at 390 MHz (forward 
ES) and 11.01 GHz (backward ES). 



netic energy in the phonon wavepacket, w 2 /2 J s p\ui\ 2 , as a 
function of frequency detuning. The agreement with exper- 
imental spectra is excellent in the backward case, with only 
one Lorentzian peak appearing. In the forward case, the ap- 
pearance of a sequence of sharp resonance peaks and its grad- 
ual disapearance with increasing frequency are correctly re- 
produced, though relative peak heights are not accurately pre- 
dicted. We display in Fig. 1(f) and Fig. 2(f) the kinetic en- 
ergy densities at resonance for selected frequencies in the for- 
ward ES spectrum. For SMF28, the energy distribution fills 
the whole fiber. For PCF, in contrast, the phonon wavepacket 
is confined by the holey microstructure. For both fibers, the 
maximum displacement is of the order of 150 fm and is dom- 
inantly transversal. The situation is dramatically different in 
the backward case, since the phonon wavepacket is found to 
be confined to the core for both fibers. For both fibers, the 
maximum displacement is of the order of 50 fm and is dom- 
inantly longitudinal, though transversal displacements don't 
vanish. 

Let us discuss the above results in the light of elastic nor- 
mal modes of the fiber section. Such an approach has been 
successful to describe GAWBS in classical fibers [16] and 
recently in PCF [2, 8, 17, 18]. It has also been used to de- 
scribe the phononic crystal properties of PCF [7, 11]. For 
a given wavevector fc, modes are obtained in the absence 
of any applied force as solutions of the eigenvalue problem 
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FIG. 3. Normalized X-profiles of optical power density and kinetic 
energy density in (a) photonic crystal fiber and (b) standard fiber. The 
black line shows the fundamental optical mode guided by total inter- 
nal reflection. The dashed blue line shows the most efficient elas- 
tic mode in the forward ES case and the red line shows the phonon 
wavepacket excited in the backward ES case. 



KU = cu 2 MU. Because the fiber cross-section has finite di- 
mensions, the spectrum of eigenmodes u) n is discrete. Modes 
are such that their potential and kinetic energy are equal, 
U^KUn — ui^U^MUn. They are furthermore orthogonal, 
or U} n KU n — U^MUn = 0, provided u m ^ oj n . They 
thus constitute a complete basis, and the solution to the forced 
equation can be expanded as 



U ■ 



E 



U^KU n 



(6) 



As a consequence, the phonon response to the electrostric- 
tive excitation is a sequence of Lorentzian-shaped reso- 
nances; at each resonance frequency, the modal shape of the 
phonon wavepacket is a normal mode of the fiber. Clearly, 
the electrostriction gain is governed by the overlap integral 
UlXT^ = J v Ar(u n )* d T^, that compares the strain field 
of each mode with the electrostrictive stress distribution. Be- 
cause of elastic losses, the response cannot diverge to infinity 
at resonance. Writing u) n = us rn + iu)i n , 



U(lo 



(7) 



The first factor in this expression is proportional to the Ce- 
faclor of the resonance. As the Q-factor in our viscoelastic 
loss model scales with the inverse of frequency, it is clear that 
resonance peaks must decrease with increasing frequency, as 
observed experimentally and with the direct computation in 
the forward ES case. 

Now it is clear that the above normal mode argument cannot 
explain the observations for backward ES (fc = 2ki). Actu- 
ally, as k becomes large, there exist more and more normal 
modes within any frequency interval, and the above analysis 
in terms of a discrete spectrum of modes becomes inadequate. 
The backward case can instead be treated under the following 
short-wavelength approximation. For an optical wavelength 
of 1.55 /im in a vacuum, indeed, the optical wavelength in sil- 
ica is about 1 /im and the acoustic phonon wavelength is thus 
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about 500 nm. This is typically much smaller than the diame- 
ter of a classical optical fiber, but still comparable to the size 
of the core of some PCF or integrated optical waveguides. We 
consider specifically core diameters exceeding a few microns 
in the following discussion. Neglecting Kq and kK\ in the 
stiffness matrix expression, we have 

k 2 K 2 U B =uj 2 b MU b , (8) 

which defines the usual Brillouin frequency wg. Ub, the 
eigenvector corresponding to eigenfrequency ljb, is a sim- 
ple plane wave [19]. In the short-wavelength approximation, 
X ps kXi and the electrostrictive force is mostly longitudinal 
because of the particular symmetries of the elasto-optic tensor 
of silica. As a consequence, shear waves satisfying Eq. (8) 
will hardly be excited and we can take Ub as dominantly lon- 
gitudinally polarized, and the velocity ug/is as very close to 
the velocity for longitudinal waves in silica, \fc\xfp = 5950 
m/s. It must be stressed that this situation is peculiar to silica, 
and would not necessarily apply to other materials, such as 
silicon. We further observe that because silica is isotropic, the 
solution to Eq. (8) is actually independent of the propagation 
direction. The actual response of the medium to the ES stress, 
U, can then be decomposed as an angular spectrum of plane 
waves, satisfying k 2 KiU = ui B MU to a good approximation 
for to ps u>b- Replacing this relation in Eq. (5), we obtain that 

U ps [K + k B K x + (w| - u^M^XT*. (9) 

This expression represents a Lorentzian distribution centered 
on the Brillouin frequency. Close to resonance, the term 
(lu 2 b —oj 2 ) is minimum, but non vanishing because of phonon 
loss. Significantly, the transverse distribution of the phonon 
wavepacket is similar to the electrostrictive stress distribution, 
T es . This is very different from the forward case, for which 
this distribution is given by the normal modes of the fiber, and 
thus extends throughout the fiber cross-section. Fig. 3 illus- 
trates this point by showing profiles of the distribution of the 
phonon wavepackets. We believe this particular distribution 
of the phonon wavepacket explains why SBS is remarkably 
efficient in optical fibers, independently of the fact that the in- 
volved acoustic phonons are guided or not by the transversal 
structure of the fiber. 

As a conclusion, the consideration of electrostrictive forces 
has allowed us to evaluate precisely the phonon wavepackets 
generated in optical fibers by Brillouin-type interactions. Our 
model encompasses the forward interaction case, including 
GAWBS, and the backward interaction case, including SBS. 
In all expressions we have derived, the frequency response to 
electrostriction is a Lorentzian function whose width is dic- 
tated by intrinsic phonon losses, or the inverse of the phonon 
life time, possibly dependent on frequency. Significantly, the 
phonon wavepacket generated via backward ES is naturally 
guided by the light that gave it birth. Our method could fur- 
thermore be extended to almost any type of optical waveg- 
uide. Silicon nanostructures, for instance, have attracted at- 



tention recently in view of obtaining optomechanical cou- 
pling in structures supporting simultaneous confinement of 
elastic and optical waves [20]. Optical forces, including ES 
and radiation pressure [21], have been predicted to scale to 
large values in nanoscale waveguides [3]. The model we have 
presented, possibly supplemented with optical surface force 
terms, should yield the phonon wavepackets generated in such 
cases. 
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